import sympy
from sympy import *
from sympy.abc import x,y,z,f
import numpy as np
f = (3*x**2+y**2)**(4*x+2*y)
fx = diff(f,x)
fy = diff(f,y)
print(fx)
print(fx.evalf(subs={x:1,y:2}))
print(fy.evalf(subs={x:1,y:2}))